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Abstract. We consider spherically symmetric inhomogeneous dust models with 
a positive cosmological constant. A, given by the Lemaitre-Tolman-Bondi metric. 
These configurations provide a simple but useful generalization of the A— CDM 
model describing cold dark matter (CDM) and a A term, which seems to fit current 
cosmological observations. The dynamics of these models can be fully described 
by scalar evolution equations that can be given in the form of a proper dynamical 
system associated with a 4— dimensional phase space whose critical points and 
invariant subspaces are examined and classified. The phase space evolution 
of various configurations is studied in detail by means of two 2-dimensional 
subspaces; a projection into the invariant homogeneous subspace associated with 
A-CDM solutions with FLRW metric, and a projection into a subspace generated 
by suitably defined fiuctuations that convey the effects of inhomogeneity. We 
look at cases with perpetual expansion, bouncing and loitering behavior, as well 
as configurations with "mixed" kinematic patters, such as a collapsing region in 
an expanding background. In all cases, phase space trajectories emerge from and 
converge to stable past and future attractors in a qualitatively analogous way as 
in the case of the FLRW limit. However, we can identify in both projections 
of the phase space various qualitative features absent in the FLRW limit that 
can be useful in the construction of toy models of astrophysical and cosmological 
inhomogeneities. 



PACS numbers: 98.80.-k, 04.20.-q, 95.36. +x, 95.35.+d 

1. Introduction. 

From a phenomenological and empiric point of view, cosmological observations are 
usually tested as a first attempt by fitting them to the so-called A-CDM model, 
which is a FLRW spacetime with flat spacelike sections, whose source is dust (cold 
dark matter CDM) and a A field (dark energy) (see pQ |2] for comprehensive reviews) . 
Although the universe may be nearly homogeneous at scales larger than the so-called 
homogeneity scale (over 150-300 Mpc), thus justifying the use of linear perturbations 
in its dynamical study, it is clearly inhomogeneous at smaller scales in which structure 
formation mostly involving CDM has taken place. Hence, the study of inhomogeneous 
sources made of dust and a nonzero "A field" is a very relevant topic. In particular, 
spacetimes of this type of source with spherical symmetry provide simple but non- 
trivial inhomogeneous generalizations of the A-CDM model. 
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Inhomogeneous spherically symmetric dust solutions of Einstein's equations with 
A = 0, described by the Lemaitre-Tolman-Bondi (LTB) metric have been widely 
used to construct simple models of cosmological inhomogeneities (see [31 3] for a 
comprehensive discussion and review). However, the literature on LTB models with 
A > (which we will denote by A-LTB models) is much less abundant E] , dealing 
mostly on singularities and censorship [7] . See [8] for a recent study of the dynamics 
of LTB models (considering also A 7^ 0), while Szekeres models with A ^ have been 
examined in 

In the present paper we attempt to study the A~LTB models under the qualitative 
and numerical phase space techniques known as "dynamical systems" . These 
techniques have been applied successfully to a wide variety of spacetimes (FLRW, 
Kantowski-Sachs, Bianchi models, self-similar and static spacetimes, etc), for which 
evolution equations for the "expansion normalized" phase variables can be reduced 
to a system of autonomous ODE's that define a self-consistent phase space (see 
[TOl nn IT^ IT5] for a comprehensive elaboration and review) . A dynamical system 
analysis has been been successfully applied to LTB dust solutions with A = in 
reference [M]. We generalize here this study to the case A > 0. The plan of this 
article is summarized below. 

The class of A-LTB spacetimes and their corresponding covariant objects are 
introduced in section 2, while in section 3 we show that their dynamics can be 
completely determined by "quasi-local" (QL) variables [13 (Ml HZj that lead to a 
set of "fluid flow" evolution equations |18j and an initial value parametrization that 
is useful for carrying on numeric work on the models. These evolution equations 
are transformed in section 4 into a 4-dimensional dynamical system constructed by 
"expansion normalized" variables [TUllTTlfT^llSj . In section 5 we list the critical points 
(past and future attractors and saddle points) of this dynamical system and classify 
all the invariant subspaces of the phase space: the homogeneous FLRW subset, the 
spatially flat subspace, as well as the Kottler vacuum solution 19, (Schwarzschild-de 
Sitter) and "pure" dust with A = subspaces. We examine in section 6 the different 
patterns of kinematic evolution of A-LTB models. The admissible topological classes 
of the space slices orthogonal to the 4-velocity are summarized in section 7. Since 
the phase space is 4-dimensional, we decompose it in section 8 in two 2-dimensional 
projections: the "homogeneous" and "inhomogeneous" subsystems. These projections 
convey the full dynamical information, and thus are used in section 9 to study in detail 
the phase space trajectories of representative configurations associated with each one 
of the kinematic patterns and topologies of the space slices listed in sections 6 and 
7. For this purpose, we solve numerically the dynamical system obtained in section 4 
for each configuration, together with the QL evolution equations derived in section 3 
(which is necessary when the phase space variables diverge in collapsing or bouncing 
configurations). Section 10 summarizes the results and provides a final discussion. 

2. A LTB spacetimes 

Spherically symmetric dust solutions with nonzero cosmological constant. A, or (as 
we shall call them) "A-Lemaitre-Tolman-Bondi" (A-LTB) models are characterized 
by the LTB metric [3 Sj 

ds^ = -c^dt^ + ^^-^ + i?2 [de^ + sin2 9 d<p'] , (1) 
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where R — R{x°, r), (a;° — ct), R' — dR/dr, E = E{r), and the energy-momentum 
tensor 

^"'' = pM"M''-A5°^ (2) 

where m° = Sq and p{x''\r) is the rest-mass density. The field equations G"^ — kT"^^ 
(with = SttG/c^) for yield: 

KQpc^R^R' = 2Af', (3) 



i?^ = ^ + f Ai?^ + E, (4) 

where where R = dR/dx^, and the function M = M{r), which appears as an 
"integration constant', is the "effective" rest mass-energy in length units. 

The comoving geodesic 4-velocity field defines a foliation of spacelike 
hypersurfaces, ^T{t), orthogonal to u", marked by t constant, and having an induced 
metric hat = "aUb + 9ab — ^a^igij^ whcrc i,j — r,9,4>. Besides p and A, the 
remaining covariant objects associated with A-LTB spacetimes are two scalars: the 
expansion, Q — VqU°, and ^TZ^ the Ricci scalar of hypersurfaces ^T{t): 

i? R^R' ' 

plus two spacelike traceless tensors, the shear tensor, aab — (a^b) — (9/3) hab, and 
the electric Weyl tensor, Eab — Cacbdu'^u'^- Both of these tensors can be completely 
and covariantly determined in terms of a single scalar function as fx"'' = S S"'' and 
^ab ^ g^ab^ where S"^ = h"'' - Sx'^x'' with x" = VhT^ 6^ being the unit vector 
orthogonal to and to the orbits of S0(3). The functions E and £ are 
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R R' 

r'^ W 



f = -ypc+^. (6) 



The set {p. A, 6, S, £, ^TZ] then provides a complete representation of local covariant 
scalars for the A-LTB models. 



3. Quasi local (QL) scalar variables. 



Following dUdS], we define for the scalars A = {p, 9, ^TZ, A} their dual quasi — local 
QL scalars Aq — {pq, 9,, ^TZq, A,} by means of the integralm 

J^^I^QAR^R'dx 

iZoR'R'dx 

where the integral is evaluated in spherical comoving domains P = x i?[r] C ^T(t), 
where is the unit 2-sphere parametrized by (9, 4>) and i?[r] = {a; | < a; < r}, where 



-^q cx—r ^9 J • ('^) 



f We are using a different notation from that of I15II16| . where QL functions Aq were denoted by A». 
These same functions where introduced in [14| . but were not called QL functions. In that reference 
the notation used was: £, H, S respectively for L, 69/3, -S'-"') and A'™), A^*) for S^^*), while 
(p), i'^'R-) and (Q) were used for pq, ^Tiq and The use of the notation (..) can be misleading 

because |7| is not an average unless it is defined as a functional (see 1201 ). See also 1151 1161 121| 
for a discussion of the covariant nature of the integral Q and its relation to the Misner-Sharp QL 
mass-energy function. 
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X — marks a symmetry center. It is evident that A, = A. [|] From ^ and ([t]), the 
QL scalars are 



Ko 



■Pi 



A„ = A. 



2M Qq _ R ■^TZq _ E 

The local and QL scalars can be related by means of "relative deviations" 



A„ 



A'q/Aq 



(8) 



(9) 



.q 3R'/R 

where we used the property: A'^ — {3R'/R){A — Aq) which follows directly from ([t]). 

It is convenient for numerical work to use dimensionless variables, hence we 
introduce 

" 'TZq ,_^oA 

6hr 



2m< 



_ KoPqC 



A = 



hs c(t - <i), 



(10) 



where hj^ is an arbitrary length scale and t = t^ marks a fiducial initial hypersurface 
^71, so that it coincides with t — 0. Considering the definitions ( 10 ) and following 
[T^ [TB] , the dynamics of LTB-dS models is fully determined by the following evolution 
equations [U 



djiq 
dr 

dZq 

Ih 

dr ' 
dr ^ 

together with the constraint 



3^q : 



'9' 

P-q + A, 

1 + J(^) 

3<5( 



1 



^iq -A(5(^) 



2/ig — Kg + A, 



Once the system (lla)-(ll(i) is solved, all scalars {p, A, 6, S, £, ^7?,} follow as: 



1 + 5^^) 



1 + 5^^) 



l + J^'^' 



(11a) 
(116) 
(11c) 
(lid) 

(12) 

(13a) 
(136) 



i: = -hsZq5^^\ 

As proven in |15j . evolution equations like (lla|-(lld| completely determine the 
dynamics of all spacetimes described by an LTB metric in the comoving frame 
(including A-LTB models). In this description the variables Aq (which satisfy FLRW 
dynamics) formally define a FLRW background, while the (5^^^ can be rigorously 
characterized as covariant, gauge-invariant and non-linear perturbations on this 
background. 

It is useful to introduce an initial value formulation for the LTB-dS models. By 
considering an arbitrary slice "^Ti — '^T[t,;], we rephrase the metric functions R and R' 
in ([T]) as dimensionless scale factors 

^, r = 4/^=i + 4/^. (14) 
Ri R'i/Ri R'i/Ri 



L = 



§ The radial range can bo modified to include models not admitting a symmetry center or when the 
^T[t] are res tricte d by curvature singularities, see [151 1161 ITt] 
II Equation (lldi was given incorrectly in |15l I16| . 
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where the subindex i will denote henceforth evaluation ait — ti. Considering as initial 
conditions for (lla)-(lld) the initial value functions {jiqi, Hqi, A}, the standard free 



parameters in ( 3 ) follow from ([s]) as M = /i^i h'^ i?f and E 
from Q, Q and (13 a I to the Friedman-like equation 

Q{L) 



-Kqi hi. Rf, leading 



L 



Q{L) = 2Ai, 



XL' 



and the following scaling laws 

l + S. 



1 + = 



(m) 



r 



where the "Omega" factors are 



^[p) ^ ^0 ^ 2^ 



1 - - 



3 /ig 
A 

3/^2^2 



2;; 



A 2XL' 

Z2 ^ 



Q{L)' 



(15) 

(16) 

(17) 
(18) 

(19) 
(20) 



The metric ([ij in this initial value formulation takes the following FLRW-like form: 
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r2 i?'2dr2 



1 



Kqi hi i?j 



SIT? ed(j)^ 



(21) 



so that L{t, r) = and T{t, r) = Q respectively mark the central and shell crossing 
singularities. The scale factor L follows from the implicit solutions of ( [T5| , which 
as opposed to the case A = 0, cannot be given (in general) in terms of elementary 
functions, but by elliptic integrals that can be formally written as |S] 

(22) 



T — Z — Zi, Z — Z (L, flqi, Kqi, A) — 



2flqi + AL3 - KqiL] 



1/2 ' 



Hence, a formal expression for T can be obtained by implicit radial derivation of ( 22 ) 

r 



1 



3[^qj {Z - Zi) 



+ Kqi (Z - Zi)^i^^. S. 



LZ_, 



(23) 



where the subindices and denote partial derivatives with respect to ^qi and 
Kqi, while Zi is the elliptic integral Z evaluated at r = (L = 1). As a consequence, 
the conditions to avoid shell crossings, F > 0, cannot be (in general) given analytically 
as in the case A = (see [331 122] ) ■ 

^ The metric ( |2l| simplifies with the coordinate choice ( |46| in section 7, leaving hl^^ as the only 
characteristic length. 
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4. A dynamical system. 



The evolution equations (lla)-(lld) lead in a natural manner to a dynamical system 
based on "expansion normalized" variables (TUl El HH [13] . Following the methodology 
of [14] , we define the dynamical system variables by normalizing with the QL expansion 



scalar, Zq, and by considering the "Omega" parameters given by (19|-(20). We define 



also a suitable evolution parameter ^ by means of the coordinate transformation 

r = r(^,f), r = f, (24) 

so that for every scalar A(t, r) = A{T{^,r),r) — A{^,r). Choosing ^(r, r) = ln(i), 



so that [d^/dT]r — Zq 

differential operator 

d_ _ 



hg ^R/R, leads for comoving curves 



_3_d_ 
Qqcdt' 



constant) to the 



(25) 



Eliminating /i^ and A in terms of rj''^' and Cl^^^ from (19)-(20l and using (25 1, the 



system (lla)-(lld) transforms into the following dynamical system 



= -3 



1 + ^(2) 



1 + 3(5(^' 



(26a) 
(266) 
(26c) 

Cl^^U^'-\ (26d) 



which can be integrated by specifying initial conditions at a suitable initial 
hypersurface marked by constant ^. Notice that the hypersurfaces of constant ^ do 
not coincide with the marked by constant t 7^ (though the initial hypersurface 
r = coincides with ^ = 0). Hence, the radial gradient of a scalar A along slices with 
constant t 7^ are different from those along slices with constant ^ 7^ 0: 



A' = 



dA' 




'dA' 




'dA' 




'dA' 


dr 


T 


dr 


^ dr 




r 


dr 



dA 



dr 



(27) 



where we usedj 25 1 and ( 23 1 to eliminate ^' and 9/9^ in terms of F and d/dr. Applying 
(271 to (19)-(20) leads to the gradients of and Jl'^-* along surfaces ^ = constant 









dr 




Hi 












dr 




Ri 



F - 1 
3F 



CiM _ 2fi(^) - 1 



$7(^)F 



F - 1 
3F 



(28a) 



(286) 



which provide consistency conditions for (26a)-(26d), as the mixed derivatives d'^/d£,dr 
obtained from these equations coincide with those obtained from the system. 



Equations (26a\-{26d\ are autonomous evolution PDE's that contain only 
derivatives with respect to ^, and so the radial dependence of the functions enters as 
a parameter. Still, we are dealing with PDE's, while the usual dynamical systems 
techniques (as discussed in [TOl [lU [121 [13]) bave been conceived for autonomous 
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Table 1. Critical points of various invariant subsets. The letters PA and FA 
stand for past and future attractors (source/sink), while SP is a saddle point. See 
l |31| |-| |35l l for reference of these critical points (except those of the Homogeneous 
subspace). 



Subspace 


Critical 


points 












Homogeneous 




= 1, 


= 


0, 


(5(^) = 




= 0, 


hPA 


subspace (FLRW) 




= 0, 


= 


1, 


(5(^) = 




= 0, 


hFA 


= = 




= 0, 


= 


0, 


(5(^) = 




0, 


hSP2 


Spatially flat 




= 1, 


o(^) = 


0, 


(5(^) = -1, 




--1/2, 


PA 


subspace 




= 0, 




1, 


J^'*-' arbitrary, 


S(-) 


= 0, 


FA 


+ nw = 1 




= 1, 


f](^) = 


0, 


= -1, 


S(-) 


= 1/3, 


SP3 






= 1, 


o(^) = 


0, 


(5(^) = 0, 


Si-) 


= 0, 


SP4 






= 0, 




1, 


(5'^) = -1, 


Si-) 


- -2/3, 


SP5 


Schwarzschild- 




= 1, 


= 


0, 


(5(^) = -1, 


Si-) 


= -1/2, 


PA 


deSitter subspace 




= 0, 


iiw = 


1, 


(5(^) = -1, 


Si-) 


= 0, 


FA 


J^'^) = -1 




= 0, 


iiw = 


0, 


(5(^) = -1, 


Si-) 


= -1/3, 


SPl 






= 0, 




0, 


(5(^) = -1, 


Si-) 


= 0, 


SP2 






= 1, 


f](^) = 


0, 


(5(^) = -1, 


Si-) 


= 1/3, 


SP3 






= 0, 




1, 


(5(^) = -1, 


Si-) 


= -2/3, 


SP5 


LTB dust 




1, 


f2(^) = 


0, 


(5(^) = -1, 


Si-) 


= -1/2, 


PA 


without A 




= 0, 


f2(^) - 


0, 


J^'^'' arbitrary, 


Si-) 


= 0, 


FA 


A(^) = 




= 0, 


= 


0, 


(5(^) = -1, 


Si-) 


= -1/3, 


SPl 






= 1, 


= 


0, 


(5(^) = -1, 


Si-) 


= 1/3, 


SP3 






= 1, 


ii(^) = 


0, 


(5(^) = 0, 


Si-) 


= , 


SP4 



systems of ODE's. However, as shown rigorously in Appendix B of 14 , a system 
of PDE's like (26a)-(26d) can be considered as a proper dynamical system, but one 
that is constrained by the fulfillment of the the initial conditions given explicitly at 
^ = 



(m) 



Ci':'^ = 




(29a) 
(29&) 

(29c) 
(29 d) 



where we used ^ and (19)-(20l evaluated at ^ = (so that L = F = 1). Once these 
conditions are satisfied in the initial hypersurface f = (which coincides with t = 0), 
then we can apply to the numerical solutions of ( 26a[ )-(26(i) the standard techniques 
of dynamical systems. 
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5. The phase space, invariant subsets and critical points. 

5.1. The phase space in general. 

The solution curves [fl'^) (C, r), r), r), (5'^) (C, r)] of the dynamical 



system (26a)-(26(i|, obtained by numerical integration for initial conditions (29a 



(29 d\, evolve in a 4-dimensional phase space parametrized as 

V = Cl'-^^ S^^'\ c R^. 

and characterized by the following critical points: 
• Source (past attractor): PA 

j(^) =-1/2, 



-1, 



n^'^) = 1, n^^'> = 0, s^''^ === 

• Sink (future attractor): FA 

(5^'^-' arbitrary, 



• Saddle Points: SP 



5<^^'^ = -1, 
(5('') arbitrary. 



1, n^^^ = 0, 



(5('^) = 0, 

fi(^) = 0, n^^'> = 1, = -1, 



<5(") = 0, 

S(^) = -1/3, 
= 0, 

<5(^) = 1/3, 
<5(^) = 0, 

= -2/3. 



SPl 
SP2 



SP3 
SP4 

SP5. 



(30) 

(31) 
(32) 

(33) 

(34) 
(35) 



As we show in the following subsections, these critical points are common to various 
invariant subspaces. Notice that the future attractor is a set of sinks defining a curve 
in V parametrized by 5'^''^ (the same remark applies to the saddle and the saddle 
SP2). 



5.2. The homogeneous subspace and symmetry centers: FLRW dynamics. 



Bearing in mind (17), (|T8|) and (p3|, we have: 



which defines the "homogeneous" invariant subspace: 
H EE 0, 0} C V, 



V e 7^ 0, 



(36) 



(37) 



characterized by the 2-dimensional dynamical system (26a)-(26(i), which only involves 
rjf'^' and tl'^^^ , and is formally identical to the dynamical system associated with FLRW 
spacetimes with a dust plus A source. Analytic solutions for (26a)-( 26(j| are given by 
(191 and (20) with L — exp(^), and the critical points are listed in table [T] they are 
the projection of PA, FA and SP2 on (37). Notice that, in general and irrespectively 



of initial conditions, the following limits hold for all curves in H: 

C ^ oo (or L oo) then fl^^^ 0, O!-^^ 1, (38a) 

C -oo (or L ^ 0) then f^^'"' ^ 1, VL^^^ ^ 0, (386) 

Hence, the past and future attractors of H clearly correspond to the past and future 
limits ^ — >■ ±cx) of the inextensible trajectories in given by (19)-(20). 

If r = marks a symmetry center, then (36) must hold for all ^ at r = 0. Hence, 



irrespectively of the initial conditions for r > 0, if a center worldline exists it is the 



only trajectory in V whose phase space evolution takes place entirely in (37). This is 
also valid if there is a second symmetry center marked by r = Tc > 0. 
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5.3. The spatially flat sub space. 



Bearing in mind the forms of r^'^^ and $7'^'^-' given by ( 19 )-( 20 1, we obtain the fohowing 
important resuh 

+ f2f ^ - 1 = ^ n^"^ + n^^'> -1^0, v^T^o. (39) 

Since Cl'^'^^ + Cl'-'^^ = 1 is equivalent to Hq = 0, and Kq = implies '^TZq = E = 0, which 



(from (|5|) implies ^TZ = 0, the constraint (391 defines the spatially flat invariant subset 



SF = f2(^\ (5^^), (5(^) I f2(^) + fi(^) - 1 = 0} c 7^, (40) 

Notice that the intersection of SF and H is the straight line {ri(^)(^) + Cl^^\£,) = 1} 
in the [f]*^^-*, 51'^'''^] plane. This is the invariant subspace that corresponds to the "A- 
CDM" model: a FLRW cosmology whose source is dust (CDM) and a cosmological 
constant. For Jl'^^-' ^ 0.3 and Cl^^'> ^ 0.7, this model is often regarded as the best fit 
to cosmological observations at large scales and in the present cosmic era. 

The 3-dimensional dynamical system associated with SF follows from the general 
system by the restriction Sl'^'^^ = 1 — Cl^^'' . Its critical points are listed in table [l] 
Notice that the past and future attractors exactly coincide with those of the general 
case. This indicates that all models have an initial asymptotic state with zero spatial 
curvature, while all models in which ^ can be infinitely extended evolve towards a final 
zero curvature asymptotic state. The saddle points SPl and SP2 are not critical 
points. 

5.4. The Kottler solution: vacuum Schwarzschild-deSitter subspace 
Consider initial conditions with arbitrary tl\^\d^'^ but (J-^' = —1. Then, if 



there are no shell crossings (F > 0) the scaling law (HTI implies that (5'^'(0 = -1 
holds for all ^. By inserting (5^'^' = — 1 in ([s]) and considering ([3|, this constraint 
defines an invariant 3-dimensional subspace 

V = {f7(^), f2(^\ -1, (5(^)} CT', (41) 

characterized by: p = 0, with {hqc^ /i)pqR^ = 2^qhlR^ = 2Mq and /Lt, = Mq/R^. 
These parameters identify the vacuum Kottler solution [19] (or Schwarzschild-de 
Sitter) as the vacuum limit of A-LTB models with Schwarzschild mass Mq, and given 
in comoving coordinates constructed with radial geodesic observers (notice that Q 
with M — Mq is the equation of these geodesies). 



The 3-dimensional dynamical system in this case follows by substituting ( 41 ) 
into (26a)-(26(i), while the corresponding critical points are listed in Table [T] Notice 



that the source PA and all saddles except SP4 exactly coincide with those of the 
general case, whereas the future attractor FA is contained in the future attractor of 
the general case: it is the point (5'^' = — 1 in the curve of sinks parametrized by S^'^^ 

5.5. LTB dust without A. 

If initial conditions are selected with A = = f)f' = 0, then Vt^^\0 = holds for all 
^. In this case we have the invariant 3-dimensional subspace 

D = 0, J^'^), C (42) 

associated with the standard dust LTB solutions without cosmological constant. This 
subset was examined in detail in |14) The corresponding 3-dimensional dynamical 

+ The phase space variables fj'''', 5^^^"> and (5'^^ were denoted in [H] by {Q), A(™) and -S. 
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system follows by setting fi^^^ = in (26a)-{26d), with the critical points listed in 



Table [TJ Notice that the source exactly coincides with the source PA of the general 
case, indicating how the cosmological constant plays no role in the early evolution 
(near the initial curvature singularity ^ —oo). The saddle points SP2 and SP5 are 
not critical points of this subsystem. The future attractor FA is different from that of 
the general case, which indicates that A plays a distinctive dominant role in the late 
evolution in layers in which ^ — > oo. 



6. Kinematics of dust layers. 



The relation between the possible type of kinematic evolution and initial conditions 



can be appreciated by looking at the roots of L^^ in equation (151, or equivalently, by 



looking at the positive real roots of the cubic Q{L), which we rewrite as 
Q{L) ^ P.L'' -a,L + l = Q, 

Cl^^^ A i]^^^ + f2^^) - 1 



ft 



(m) 



2^q 



(43 a) 
(436) 



While (43 a I can be solved exactly, it is more practical to examine it by looking at the 



qualitative behavior of the cubic polynomial Q{L) (see figure [Tp). The following cases 
arise 



Q(L) 




(b) 



X - dominated 



Qmiii<0 X,^ 





dominated 



<0 Q(m) 



Figure 1. Initial conditions and Kinematics. Panel (a) shows the 
cubic polynomial Q(L) defined in ( |43a[ l, whose roots determine the zeroes 
of Lt-. Configurations corresponding to Qmi„ > (between the two thick 
curves) are perpetually expanding (no zeroes) but have positive spatial curvature. 
Configurations with Q„ii„ < must bounce or collapse. Panel (b) shows the 
regions in the space of initial conditions {O^'*' , A'^'} for which Qmi„ < holds 
(shaded areas). 



If ai < then Q{L) has no extrema and Q{L) > 1 holds for all L. This case 
clearly corresponds to the perpetually expanding/collapsing configuration with 
Ci^+qW < 1 (or < 0). 
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If Q!i > then Q{L) has a minimum at L = i,„i„ = [f3i/{3ai)Y^^ , so that 



K 



3/2 2 



3/2 



Q(L_) = 1 - = 1 - ^ . . ' ■ (44) 



(A) 



(m) 



Clearly, since QiL) > then (43 a I has no positive roots if Q,„i„ > (see figure 

[T^), which corresponds to perpetually expanding/collapsing configurations with 
Cl['^^ + flp^ > 1 (or Kqi > 0). However, a sufficient condition for the existence of 



positive roots of (43 a 



IS 

Q_(f2(^\$7f') <0, (45) 

which corresponds to the configurations that have turning values in which the 
range of L is restricted by a bounce where L r ~ 0. 

In order to examine the initial conditions associated with the possible kinematic 
evolutions described above, we examine the constraint (451 by plotting in figure [TjD 



the level curves of Q,„i„ in the parameter plane {ri['^\Cl^}. Since these initial value 
functions depend on r, then any given set of initial co nditions will appear as a curve 
['tll''\r),Clf'\r)] in this plane. The regions where holds are marked by the gray 
areas in the upper left and lower right corners in figure [TJa, respectively corresponding 
to a "A-dominated" and "^-dominated" regimes. It is interes ting to notice that the 
"/x-dominated" regime is more restrictive for the fulfillment of ( 45 ) , since a larger Cl^'^'^ 

is needed even for smaU Cl[^^ < 1. 

The kinematic evolution patterns of dust layers in A-LTB models and their 
relation with initial conditions are described below: 

(i) Perpetual expansion (collapse) from (to) a singularity. If Kqi < 0, or if Kqi > 

but Cl^^'^ and Cl^^"^ are roughly comparable for all r, we are well within the non- 
shaded area of figure [I]d and ( 45 ) does not hold. See figures [2] [s] and [i] 



(ii) All layers initially expand from and then collapse to a singularity. This pattern 
corresponds to the extreme "^-dominated" regime in figure [T]d and to the range 
restriction < L < Li with the bounce L^r = at L = Li. Initial conditions 
must comply with > 1 and Cl\^^ < 1 for all r. If il^^' = (pure dust) then 
Kqi > is a sufficient condition for this evolution, but Kqi > is only necessary 
(not sufficient) if Cl'-^^ > 0. See figure [s] 

(iii) All layers initially collapse from infinity, bounce and re-expand. This pattern 
corresponds to the "A-dominated" regime in figure[l}D, and to the range restriction 
L > L2 where the bounce of L occurs at L = so that layers avoid the collapse. 
Initial conditions must comply with Kqi > 0, Cl'f'^ <^ 1 and Clf''^ ^ 1 holding for 
all layers. See figure [6j 

(iv) Loitering evolution. Layers initially expanding or collapsing reach an unstable 
static state. Layers may collapse or expand at later times (A-dust FLRW 
models with this kinematics are known as Eddington-Lemaitre models). Initial 
conditions associated with this pattern comply with Kqi > and Q„i„ = 0. See 
figure [7j 

(v) Mixed evolution patterns. Since f2|^\ f2|^^ depend on r, dust layers may 
show different evolution patterns in different radial ranges (as opposed to the 
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homogeneous FLRW case). In these cases initial conditions must be specified so 



that (451 only holds in part of the radial range, or equivalently, with the curve 
['tl\^\r) , Cl\^\r)] passing from the non-shaded area to either one of the shaded 
areas of figure [TJs. Two possible configurations are: 

• Layers around the center follow pattern (ii) and "external" layers evolve as 
in (i). See figure [Sj 

• Layers around the center follow pattern (i) and "external" layers evolve as 
in (iii). See figure [9j 

We examine the phase spaces of configurations following the evolution patterns (i)-(v) 
in section 9. See |24| for detailed plots of the scale factor for patterns (i)-(iv) in FLRW 
cosmologies. 

7. Topology of the space slices and a choice of radial coordinate. 

The LTB metric (in both forms ([!]) and (21)) admits an arbitrary rescaling of the 
radial coordinate. We use this coordinate gauge freedom to fix the function Ri in ( 14 1 



and (21 ) by the choice 

R^ = h-'f{r), (46) 



where / is a smooth dimensionless function, while hj^ (defined in (lOl) provides a 
length scale for the radial coordinate. The topological class of the slices '^T[t] is 
given by the number of symmetry centers, which determines the allowed forms of /. 
Assuming absence of shell crossings, so that F > holds everywhere, the admissible 
topologies of the slices '^T[t] are 

• "Open models" . The ^T[t] are homeomorphic to M'^ and admit only one symmetry 
center at r = 0. Since R' > and R'^ > must hold for all r, then / can be any 
monotonously increasing function complying with /(O) = 0. The simplest choice 
is /(r) = r, but we wiU use instead /(r) = tanr so that i?i — )• oo corresponds to 
the finite limit r — > n/2. 

• "Closed models" . The ^T[t] are homeomorphic to S'^ and admit two symmetry 
centers, at r = and r = rc- Since R{t, 0) = R(t, Tc) = for all t, then there must 
exist (at least) a turning value r = r^^ so that R'{ct,rt^) = where < r^^ < Tc- 
The function / must comply with /(O) = f{rc) = and /'(r^^) — 0. A convenient 
choice is /(r) = sinr. 

• "Wormholes" . The ^T[t] can be homeomorphic to x M or x §i. While 
the ^T[t] do not admit symmetry centers, there must be (at least) a turning 
value R'{ct,r-t^) = 0. The function / cannot have zeroes, but must comply with 
/'(^tv) = 0. We will use /(r) — seer in the examples of section 9. 

In order to examine the relation between the topological class of the '^T[t], the sign 
of the spatial curvature and the kinematic evolution of dust layers, we remark that, 
even if F > holds, Ri must comply with the following condition associated with the 
regularity of the grr metric coefficient in ([T]) and (21 ) |161 [T71 [25] 



\l-^,,R^hl\'/' |l-^,,/2|i/2 



^ 0. (47) 
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Hence, if there is a turning value r — r^^ such that R[{ri^) ~ 0, the denominators 
in the expressions above must have a same order zero at r = r^^. This places the 
following restriction on / and Kqi 



/'(nO = ^ 



> 0. 



(48) 



If ( 47 )-( 48 1 fail to hold we have a surface layer singularity at r = r^^ [TBI [53] . 
Therefore, (48 1 implies that all regular models with closed and wormhole topologies 

1 > 0. In the case 



must have positive spatial curvature: 



> or r2. 



(A) 



of pure dust LTB models (A = 0), the implication of (48) is that dust layers of all 
models with these topologies must bounce and collapse (as > implies that L t 
has necessarily a zero). However, the situation is different if A > 0, since dust layers 
can expand for all r when Kqi > as long as the initial value functions are selected so 
that (45) does not hold (see figure [I]). Hence, fully regular and perpetually expanding 
LTB-A models with closed or wormhole topologies are possible. 



8. Phase space dynamics in terms of homogeneous and inhomogeneous 
projections. 



Since the homogeneous invariant subspace H defined by (37) corresponds to dust-A 
FLRW cosmologies, the dynamical variables f)''^-' and fl'^^^ of H become the standard 



FLRW Omega parameters, and the dynamical system for H (which is (26a) and (26&I) 
is equivalent to the dynamical system of these FLRW cosmologies (in the expansion 
normalized variables [101 [HI [HI [13] ) . We define for the latter the 2-dimensional phase 
space 

H(o) = C M^ such that H = i(H(o)) (49) 

where i is the inclusion map i : — > M"*, i{x^,x'^) = [x^, a;^, 0, 0]. 

Since il(o) is endowed with a clear physical and geometric meaning (the phase 
space of the homogeneous limit FLRW spacetime), while {S'^'^\5^^^ convey the 
deviation from homogeneity in the dynamics (and are gauge invariant covariant non- 
linear perturbations [15'), it is useful to express the general phase space V as the 
direct product 

7' = H(o)Xl, I={5(^), J(")}cM2, (50) 

so that trajectories and critical points of V can be described and analyzed by means 
of the following projection maps applied to every p = [f](^\ J*^''^, S^^^] e V: 

Homogeneous projection. 

Hh:7'-^H(o), such that Hh(p) = [Sl^^^ fJ^^^] G H(o), (51a) 
Inhomogeneous projection. 

HjiP^I, such that Hi(p) = [(5'^), (5*^'] e I. (516) 



While the use of (51a) and (51&) solves the practical problem of dealing with the 4- 
dimensional phase space V, these projections highlight interesting dynamical features 
of V which we discuss below. 



Since the dynamical equations (26a) and (266) associated with H(o) ^re completely 
independent of 5^^\ (J*^^-* (the coordinates of I), there is a one-one equivalence map 
between every FLRW phase space trajectory in H(o) and a phase space trajectory of 
V projected into H(o) by Hh. However, the initial conditions for each trajectory in 
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H(Q) are two constants: ^-nd so each phase space trajectory corresponds to 

a single spacetime configuration. On the other hand, because of the inhomogeneity of 
A-LTB spacetimes, initial conditions are given (in general) by functions of r. Hence, 
these trajectories form a one-parameter family of curves, so that each curve (when 
projected into H(o) by IIh) is equivalent to a trajectory of a distinct FLRW spacetime. 

The critical points of Hg are simply the critical points of H (see table [T]) projected 
under IIh, and correspond to the particular case spacetimes: 

nit') = 1, = Past attractor: spatially flat FLRW, (52a) 
nit') = 0, n'^^'> = 1 Future attractor: de Sitter, (526) 
n(f) = 0, Cl'-^) = Saddle: Minkowski, (52c) 

However, each one of the critical points in V listed in section 5 is also projected by 
Hh into one of the critical points of Hq listed above: 



PA, SP3 and SP4 are projected into the past attractor (52a 



FA and SP5 are projected into the future attractor (52&) 



• SPl and SP2 are projected into the saddle (52c I 

These projections clearly suggest characterizing every critical point in V as related to 
a given limiting spacetime, such as spatially flat FLRW, de Sitter or Minkowski. The 
degeneracy in projecting the critical points of V into any one of the critical points 
(52a|-(52c) occurs when we only look at the projection of the trajectories of V into 
H(o) by Hh, but it disappears when we examine the projection of these curves into I 
by means of (51&), whose coordinates {S^^^) , S^^)} convey the effects of inhomogeneity. 
Hence, a complete picture of the evolution of the trajectories in V and their relation 
with the critical points follows when we combine the description that follows from the 
two projections (51a) and (516). 

The projections (51a) and (516) will be used in the following section to analize 
the phase space trajectories of various representative A~LTB configurations associated 
with the kinematic patterns and topologies of the space slices discussed in sections 6 
and 7. 



9. Numerical and graphical examples. 



We examine, by means of numeric solutions of (26a)-(26(i) under initial conditions 



(29a)-(29d), the phase space evolution of various representative A~LTB configurations. 
By plotting the curve ['tl^'^\r) , Cl['^\r)] in the parameter space of figure [i]d, we can 
infer for any given set of initial value functions a kinematic evolution pattern (i)- 
(iv) described in section 6. The topological class of the space slices follows from the 
choice of f{r) in (46). To avoid dealing with shell crossing singularities, we verify 
numerically for every configuratio n un der study the fulfillment of F > (notice that 
r,^ = STzgS^^) follows from and ([Hf). 

Since the phase space V is 4-dimensional, we plot trajectories (curves 
parametrized by ^ with fixed r) for any given set of initial conditions in eight figures 
composed by two panels, each one containing a 2-dimensional graph corresponding 
to: 

• Panel "(a)" (left hand side) : the congruence of curves [Cl'^^) {S,,r),n'-^) {^,r)] given 
by the homogeneous projection of V defined by ( 51a[ ). 
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Panel "(b)" (right hand side) : the congruence of curves [^^^-^ (■f , r), J*-^^ (^, r)] given 
by the inhomogeneous projection of T-" defined by (51&I. 



These graphs display the congruences of the trajectories, together with thick black 
dots marking critical points (past and future attractors, PA and FA, and saddle 
points SPl, SP2, etc). In configurations that contain dust layers that collapse or 
bounce (subsections 9.2 and 9.3, figures 6-9) the phase space variables 17^'^^ Cl'-^^ 



and (5^^^ diverge as 



0. For these cases we plot the arctan of the variables, so 



that zbinfinity is brought to ±n/2. However, the trajectories cannot be continued by 



numerical i ntegr at ion o f {26a)-{26d) beyond these points. In these cases we must use 



the system (lla)-(ll(i) with the variable s = ZqS^^^ (which is bounded as Zq — 0) to 



examine the behavior of fJ^'*'^ and 6^^^ in the full evolution range. 



9.1. Perpetually expanding configurations: kinematic pattern (i). 

If initial conditions are chosen so that the curve ['tl[^'' (r) , Cl'f'^ (r)] lies entirely in the 
non-shaded region of figure [Tj3 and Zqi > for all r, the resulting configurations 
expand or collapse monotonously for all ^. Since a collapsing evolution would follow 
from an expanding one by a trivial sign inversion of ^, we will only consider the latter 
in the examples below. 



9.1.1. Negative spacial curvature. We consider /(r) — tanr in (461 in the range 
< r < 7r/2, so that i?j ^ oo as r — 7r/2. We consider first initial value functions 
given by 

l^qi = mio + ^3^^^ , (mio = 0, mil = 15.3), (o3a) 

nq^^kw+^^^^—^ (fcio = -1.2, fell = -0.1), (536) 
1 + /4(r) 

A = 0.1. (53c) 

The phase space trajectories are displayed by figure [2] which highlights the past 
attractor PA from which all trajectories emerge (and can be associated with the 
initial singularity). These curves approach the saddle point SP2 and terminate in a 
point of the line of the future attractors FA (thick gray line in panel (b)). because of 
the degeneracy in the ^''^^ coordinate in SP2 and FA, the curves in panel (b) cross 
the (5*^^^ axis at different values of (5''^^. Since ^qi — >■ as i?^ — oo, this configuration 
must be radially asymptotic to a Schwarzschild de-Sitter vacuum. This is consistent 
with the fact that in the limit r 7r/2 (which implies Ri — > oo) we have near vacuum 
conditions fi*-^' w and w (see [17]) I and thus the trajectories in panel (a) tend 
to the lower left corner saddle in the intermediary stage of the evolution between PA 
and FA. Also, the curves in panel (b) in the limit r — >■ 7r/2 reach FA with 5^^^ -> — 1, 
which characterizes the invariant subspace of a Schwarzschild de-Sitter vacuum (see 
section 5.4). 

If we consider initial conditions as in (53a)-(53c), but with mio = 0.1, fcio = 



—0.75, so that ^qi — > 0.1 and Kqi — > —0.75 as Ri — )■ cxd, we obtain a configuration 
with the same time evolution but radially asymptotic to a A-dust FLRW model 
with negative curvature [T71. The phase space trajectories (displayed by figure |3]) 
are qualitatively analogous to those of figure [2] but the different radial asymptotic 
conditions imply that the trajectories as r — ^ 7r/2 {Ri — )■ oo) do not reach as far as 
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Figure 2. Expanding configuration with negative curvature and 
asymptotic to Schwarzschild de Sitter. These phase space trajectories follow 
from initial conditions |[53al-(|53cl. 




PA 



Figure 3. Expanding configuration with negative curvature and 
asymptotic to FLRW. Initial conditions and figure conventions are as those 
used in figure [2| but with fj,qi tending asymptotically to a constant. 



the lower left corner in panel (a) and do not reach as far as (5''^-' — >■ — 1 in panel (b) 
(as expected if the model is radially asymptotic to a A-dust FLRW model) . 

9.1.2. Positive spatial curvature with open topology. We choose /(r) = tanr in the 
range < r < 7r/2 for the following initial value functions 

Mgi = "T-io + ^^^] j9™N° 1 (mio = 0.8, mil = 30.0), (54a) 
1 + }\r) 



fell — k 



Kqi ^ fcio + , (fcio = 0.7, fell = 15.5), (546) 

A = 0.75. (54c) 
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These configurations liave no equivalent in pure dust models rj'^-* = A = 0, as in 
this case positive curvature necessarily implies reversal of expansion and collapse (see 
section 6). The phase space trajectories are displayed in figure [4] As in the previous 
cases the curves emerge from the past attractor PA (big bang) , expand monotonously 
and end in the future attractor FA, but now the trajectories in panel (a) are located 
above the line tl^^^ + il'^^-' = 1 characteristic of zero spatial curvature. The curves 
now approach the saddle point SP4, whose projection coincides with the projection 
of PA in the {f2(^),r2(^)} plane (panel (a)). 

(a) (b) 




Figure 4. Expanding open model with positive curvature and open 
topology. Initial conditions are given by ( |54a| l- |54c| . See the text for further 
explanation. 



9.1.3. Positive spatial curvature with closed or wormhole topology. As in the previous 
case, these configurations have no equivalent in the pure dust case. 

• Closed topology. We take f{r) = sinr, so that the radial range is restricted by 
< r < TT, with two symmetry centers at r = 0, tt, and r = 7r/2 marking the 
"equator" of the 3-sphere S"^ where R takes its maximal value and R'{tt/2) = 
(see section 7). 

• Wormhole. We can use /(r) — seer, so that — 7r/2 < r < tt/2 and R has no 
zeroes (no symmetry centers), while i? — >■ (X) as r — > ±7r/2. This choice yields 
space slices homeomorphic to x K, where the locus r — can be identified with 
the "throat" of the wormhole where R takes its minimal value and R'{0) = 0. 

We do not display the phase space trajectories for these cases because for all choices of 
initial conditions they are very similar to those of the open models of figure [4j though 
some distinctive features are worth remarking. In the case of closed topology the 
trajectories of the two centers in the homogeneous projection (panel (a)) are identical 
to each other and to the curve for the center (r = 0) in panel (a) of figure |4] The 
trajectories in the range 7r/2 < r < tt can be uniquely mapped to those in the range 
< r < 7r/2. The curve for r — tt/2 (the "equator" of S'^) is analogous to the curve 
close to the line fl^^^ + f2*^^^ = 1 corresponding to radial infinity in the open models 
(r — >■ 7r/2 in both panels of figure [4|. In the wormhole case the curves for r < and 
r > ("left" and "right" of the throat r — 0) are identical and can be mapped to 
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each other, while the phase space curve of the throat is analogous to that of the center 
r = in both panels of figure |4] 



.2. Configurations in which expansion is reversed or halted for all dust layers. 



These configurations follow by choosing initial conditions complying with ( 45 ) , so that 
the resulting curve [f2^^''(r), Cl'l^\r)] is fully contained either in the "^-dominated" or 
the "A — dominated" regimes (shaded areas in figure [T|d). 



arctan 




Figure 5. All layers collapse. This configuration has space slices with closed 
topology and follows from initial conditions | |55a[ |— j55c| l. See the text for further 
explanation on the phase space evolution. 



9.2.1. Collapsing models: kinematic pattern (ii). If we choose jiqi « Kqi > 1 with 
A <C Hqi,Kqi for all r, we obtain strongly /i-dominated A-LTB configurations that 
are similar to the case of pure dust A = with positive curvature in which layers 
expand from an initial curvature singularity (big bang), bounce and collapse to a 
second singularity (big crunch). Notice that, in general, neither singularity coincides 
with a simultaneous hypersurface. The space slices can have closed, open or wormhole 
topologies. An example of initial conditions leading to this type of evolution is given 
by choosing / = sin r (closed topology) together with 

TOll 



l + fHrY 

fell 
l + /2(r)' 



(mio = 0.3, 
(fcio = 0.01, 



mil 



2.0), 
fell = 1.98), 



A = 0.01. 

where Kqi complies with ( 47 ) and ( 48 1 



(55a) 

(556) 
(55c) 



Since dust layers can only evolve in the range 
< L < Li (see section 6), the initial expansion reverses into a collapsing stage at a 
"bouncing point" L = Li{r) which is in general different for each r. At this bouncing 
point we have Zq = 0, and so fi^'^^, fl^^^ and 5^^^ diverge, and as a consequence, we 
cannot find the collapsing stage of the trajectories (which lies beyond this point) with 
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the numeric solutions of the dynamical system ( |26a |-( 26d\. The full trajectories in 
this case must follow from the numeric solutions of the system ( llfl)- pld| . 



Phase trajectories are displayed in figure [5] where we have highlighted in panel 
(b) a representative trajectory marked by arrows (to facilitate the explanation). The 
homogeneous projection of the phase space evolution (panel (a)) is very simple: the 
expanding curves emerge from the past attractor PA, reach the bouncing point where 
— )■ and f]^''-*, O^^-* diverge (marked by a square), the expansion reverses into 
a collapse and the curves return following the same trajectories (but run in reverse 
direction) towards PA, which is now the a future attractor in the collapsing stage. 
In the inhomogeneous projection (panel (b), see representative curve), curves follow 
different paths in the expanding and collapsing stages. All curves initially expand 
and emerge from PA, approach the saddle SP4, and head upwards to the "bouncing 
points" marked by the upper squares where — and (5'^^ — > oo for S'--^"^ variable. In 
the collapsing stage (see representative curve) the curves start in the bouncing points 
(lower squares) with 5*^^^ — > —oo for 6^^^ variable and evolve upwards to PA, which 
is now a future attractor. The fact that expanding and collapsing trajectories are 
the same (with ^ ^ — ^) in the homogeneous projection (panel (a)), but are wholly 
different in the inhomogeneous projection (panel (b)) is an important inherent feature 
due to the inhomogeneity of the models. 

We can obtain the same type of collapsing configurations with slices having 
open and wormhole topologies by means of similar initial conditions, but with 
/ = tanr, seer. 

9.2.2. Bouncing models without collapse: kinematic pattern (Hi). Strongly A- 
dominated configurations follow by selecting /i^i <C Kqi and fiqi <C A for all r. Since all 
layers start collapsing from infinite L but bounce and avoid collapsing to a singularity, 
these configurations have no equivalent in the pure dust case (A = f](^) = 0). The 
space slices can have closed, open or wormhole topologies. An example of initial 
conditions leading to this type of evolution follows from choosing / = tanr (open 
topology) and 

fiqi = mio + ^ ™" , , , (mio = 0.01, mn = 1.5), (56a) 

Kq»-fcio + — (fcio^lO.O, fell = 30.0), (566) 

A = 31.0. (56c) 

The phase space trajectories are displayed by figure [6j Since dust layers can only 
evolve in the range L > L2 (see section 6), the initial collapse from infinity bounces 
into an expansion at L = L2{r) which is in general different for each r. As in the 
collapsing case, at this bouncing point we have Zq = 0, and so fl^^), n'-^'^ and J^^' 
diverge, hence we have plotted the arctan of these variables in figure [6] As in the 



previous case, we have used the system (lla)-(ll(i) to obtain the collapsing stage of 
the evolution. 



The curve of past attractors PA is given by ([32|, and is the same as the curve 
of future attractors for the perpetually expanding configurations of figures [2] |3] and 
|4] Every phase space trajectory starts from a point (marked by PA) in this curve of 
past attractors (thick gray line) and terminates in a different point in the same curve 
(of future attractors, marked as FA). These attractors are projected into the same 
point in the homogeneous projection (panel (a)) and into two distinct points in the 
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arctan 



I 



(a) 




Figure 6. All layers bounce. This configuration follows from initial conditions 
l5^-(l56cl. 



thick gray line in the homogeneous projection (panel (b)). In panel (a) the initially 
collapsing trajectories start at PA, evolve upwards to the bounce points as — )■ 
and r2'^\ rj^^) — >■ oo (squares) and return to the same point (now a future attractor 
FA) following the same trajectory. In panel (b) (see representative trajectory marked 
by arrows), trajectories start in a point PA of the curve of attractors (thick gray line) 
and evolve downwards to the bounce points where Zq — >■ and S^^^ — > — oo (lower 
squares) and 5*^^^ variable. In the expanding phase trajectories start at a bounce 
points — >■ and 6^^^ — > oo (upper squares) and go downward towards the point FA 
of the curve of attractors. As in the previous case (collapsing layers), we also have 
a time symmetric reflection (^ — >■ — ^) in the collapsing/expanding evolution in the 
homogeneous variables, but not in the inhomogeneous variables. 

Bouncing configurations that avoid collapse can be constructed with space slices 
having closed and wormhole topologies by means of similar initial conditions, but with 
/ = sin r, sec r and Kqi complying with ( 47 ) and ( 48 1 . Their phase space trajectories 



are qualitatively similar to those of figure [6j 

9.2.3. Loitering models: kinematic pattern (iv). Initial conditions for such models 
follow from taking Kqi > and restricting the initial value functions so that (5„i„ = 
holds. Under these restrictions A-LTB models evolve towards an asymptotic static 
state, which is very unstable because of the extreme sensitivity to initial conditions 
exactly complying with (5,„i„ ~ (see section 6 and figure [T]). As a consequence, the 
numeric integration of (26a|-(26(i| for these initial conditions is extremely susceptible 
to numerical errors. 



However, if Q,, 
following closed analytic solution: 

1/2 



the Friedman-like equation ( 15 ) has the 



L 



2 tanh^ 



li A 
2 



{t — ti) — arctanh 



1 

71 



(57) 



where if = Kqi/{Z\) and we eliminated from (44). Since all phase space variables 
can the computed analytically from (57) by means of (14)-(20), we can plot the phase 
space trajectories given a choice of initial conditions. Considering 7^ = (2+7-'*)/(l+r^) 
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and A = 0.1, wc obtain the phase space diagrams displayed by figurc[7] All trajectories 
in both panels (a) and (b) begin at the past attractor PA (big bang) and evolve towards 
the static state characterized by Zq — > 0, which implies S'-^^ — > oo. 




5 PA I I.S 

(a) (b) 

Figure 7. Loitering layers. This configuration corresponds to initial conditions 
complying with Q,„i„ = 0. 



.3. Configurations with mixed kinematic patterns: kinematic patterns (v). 



Since initial value functions depend on r, we can select initial conditions so that 
different kinematic patters occur in the same configuration. In practical terms, such 
initial conditions must comply with (45) and would yield curves [Cl\'^\r), Vlf'\r)] that 



traverse the shaded and non-shaded areas of figure [TJs. Notice that these "mixed" 
pattern configurations can be easily constructed within a single inhomogeneous model, 
but cannot be set up with any single homogeneous model. 



arctanQ^^^ 




0.5 PA 1 !■? 

, . arctanQ^'^' 

(a) 



Figure 8. Structure formation scenario. Initial conditions for this 
configuration are given by | |58a| l- |58c| . 
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Figure 9. Mixed expansion and bounce. Initial conditions are given by 
ll5^-(l59cll. 



9.3.1. Structure formation scenario. Interesting configurations that can be 
associated with "structure formation" toy models follows from initial conditions in 
which (45 ) holds (shaded /i-dominated region in figure [T|3) in a radial range containing 
a symmetry center (0 < r < rb) and fails to hold (non-shaded region in figure [TJd) 
for all r > Tft. This corresponds to an inner spherical region undergoing local collapse 
surrounded by an expanding "background" . While this structure formation scenario 
can also occur with pure dust solutions, in this case the boundary r = rf, between 
the collapsing region and the expanding background corresponds to the worldline 
where spatial curvature passes from positive to negative or zero. However, if A > 
this scenario can also occur with all layers having positive spatial curvature, as a 
perpetually expanding background with positive curvature is possible. 

An example of initial conditions leading to the structure formation toy model 
follows by choosing /(r) — tanr (open topology) with 



^^qi ^ rniQ + ———-— , (mio = 0.0, mn = 20.1), (58a) 
1 + f'^{r) 

k — k 

i^qt = kiQ+ " /° , (fcio = -4.1, fell = 35.5), (586) 
l + /4(r) 

A = 0.75. (58c) 

The phase space trajectories are displayed by figure [8] All curves start their evolution 
at the past attractor PA marking the initial singularity (big bang). However, the 
"inner" layers around the center (r = 0) and the "outer" layers r > r^ follow a 
different phase space evolution. 

Since outer dust layers for a given r > rh expand perpetually, their evolution is 
analogous to that of perpetually expanding configurations depicted by figures [2] |3] or 
|4] (depending on the spatial curvature and the asymptotics of fiqi and Kqi). In the 
particular example we are considering, these trajectories (see figure |8]) have analogous 
forms to those of these three figures. From the past attractor PA they approach 
the saddle SP2 and end in the line of future attractors FA (whose inhomogeneous 
projection is marked as a thick gray line in panel (b)). 

The inner layers, on the other hand, evolve in a similar way as the trajectories 
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of collapsing configurations displayed by figure [5] they initially expand away from the 
past attractor PA approach the saddle SP4 and go upwards to the bounce points as 
Zq — >■ 0, with 0,'-'^\ fl^-^^ and S^^^ diverging (the bounce points are marked by small 
squares in both panels). The main difference with the "pure" collapse case of figure [s] 
is that the trajectories of inner layers closer to the outer perpetually expanding layers 
reach their bouncing points (see panel (b)) as f?'-^-' — > oo, whereas inner layers closer 
to the center do so as J*-^-* — — oo. The collapsing stage of these layers is not shown 
in figure|8j but it can be readily obtained from system ( lla)-( If d): it is analogous to 



that depicted in figure [5] the collapsing trajectories in panel (a) are simply the time 
reversal of the expanding ones, whereas in panel (b) they start from a bouncing point 
towards a future attractor that exactly coincides with PA (the collapsing singularity) 
but following a different path to that of the expanding curves. 

9.3.2. Bouncing /expanding models. Another mixed kinematics configuration follows 
if dust layers around the center expand from an initial singularity, while external layers 
collapse from infinity and bounce back. In this case initial conditions need to comply 
with Vt\^'' w Cl'^l^'^ near the center passing to Cl'^^'^ ^ Cl'f'^ as r grows, so that the 
curve [Cl\^\r),h^^\r)] would be passing from the non-shaded to area of figure [i]d 
towards the A-dominated shaded region. As an example of such evolution pattern we 
use initial conditions 

Mqi = '^lo + ^T-VT^rV' (mio = 0.01, mil = 15.0), (59a) 

K,^ = klo+'i^^^^, (fcio = 10.0, fell =30.0), {59b) 

X = 10.0, (59c) 

with / = tanr, leading to the phase space diagrams of figure [9] that clearly show the 
mixed evolution patterns. Layers expanding perpetually from an initial singularity 
behave as the curves in figure |4] they start from the past attractor PA, approach 
the saddle SP4 and terminate in the future attractor PA. The bouncing curves 
are analogous to those of figure [6} they originate from the future attractor PA 
(a past attractor for these curves) and evolve towards the bounce — >■ (or 
oo, j(^) — oo with S^f^^ variable), and then expand back to PA (now a 
future attractor). We only display in figure [9] the collapsing stage of this evolution, as 
the curves in the expanding stage are qualitatively analogous to the expanding curves 
of figure [6j 

10. Conclusion. 

We have conducted a comprehensive and detailed dynamical systems study of A~LTB 
spacetimes by means of covariant quasi-local (QL) variables, generalizing a previous 
dynamical systems study [T3| on LTB dust models with A = 0. Since the phase 



space for the dynamical system (26a|-(26d) is 4-dimensional, we have examined and 



plotted phase space trajectories in two 2-dimensional subsets: the "homogeneous" 



and "inhomogeneous" subspaces, Hq and I, defined by the projections (51 a I and (51&), 
and respectively displayed in panels (a) and (b) of figures [2]-[9] The decomposition 
of the phase space into Hq and I provides an important theoretical connection 
between the dynamical system study of A-LTB spacetimes and the perturbation 
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formalism presented in [151 [16]. Iii this formalism, the fluctuations {(5^^\ (5^^-'} (which 
parametrize I) are rigorously characterized as covariant and gauge invariant non-linear 
perturbations on a FLRW background defined by quasi-local scalars, which satisfy 
FLRW dynamics, and hence corresponding to the projected homogeneous subspace 
Ho parametrized by il'-^'>}. 

In LTB dust models with A = all configurations in which the space slices ^T[t] 
have closed (S'^) or wormhole (S"^ x M or S'^ x topologies, the regularity conditions 
(47|-(48) imply that spatial curvature must be positive and thus, perpetual expansion 
is not possible [IHl [HI US]- However, this restriction is no longer true when A > 0. 
While ( 47 )-( 48 ) still requires positive spatial curvature for configurations with these 
topologies, positive spatial curvature no longer forbids perpetual expansion (see section 
7). Hence, perpetually expanding and fully regular A-LTB models with closed or 
wormhole topologies are possible. 

The homogeneous projections in panels (a) of figures [2||9] show that phase space 
trajectories are qualitative analogous to those of dust-A FLRW models. However, 
each FLRW spacetime configuration would correspond to one (and only one) phase 
space trajectory in a {f]^''^ £7''*'^} diagram (because each spacetime configuration is 
uniquely determined by initial conditions given by constant As a contrast, 

since initial conditions for any single A-LTB configuration depend on r, all displayed 
curves in the diagrams in the panels (a) of each one of these figures corresponds to a 



single spacetime configuration. In other words, in terms of the projection (51a) each 



A-LTB configuration is formally equivalent to a superposed one-parameter family of 
dust-A FLRW models. 

In all A-LTB configurations (except those like the one in figure [6| the phase 
space trajectories initiate in the same source or past attractor PA, defined in (31) 
and associated with an initial singularity. This critical point can be projected into the 
past attractor of the homogeneous subspace, and can be identified with a spatially 
flat FLRW dust cosmology {n(f'> = 1, n^^^ = 0). Since J^''' = J*^' = -1, this critical 
point exactly coincides with the past attractors of LTB dust solutions with A = 
and spatially flat models with A > 0, thus indicating that the effects of A and of 
spatial curvature are negligible in the early stages of the evolution near the big bang. 
The past attractor PA can also be identifled with self-similar conditions [2], which 
is consistent with the fact that both spatially flat FLRW and LTB dust spacetimes 
with A = are compatible with a homothetic Killing vector. The past attractor (31 ) 
becomes the future attractor or sink for phase space trajectories of layers that collapse 
to a curvature singularity (big crunch). This is so, either for a collapse that is a trivial 
time reversal of a perpetual expansion (figures [2]-[4]) , or for configurations where all 
layers expand and collapse (figure [5]), or with mixed dynamics in which some layers 
expand perpetually and some collapse (figure [s]). 

Phase space trajectories of perpetually expanding layers {zg > 0) terminate in the 
same future attractor (sink) FA, defined in ([32]): figures [2]-[4] and in the perpetually 
expanding trajectories of figures [6] [8]and[9j This critical point can be projected into 
the future attractor of the homogeneous subspace (fl^^) = 0, Cl'-^'^ = 1), and is also 
the future attractor for spatially flat configurations. Hence, it can be identified with 
de Sitter spacetime and spatially fiat conditions, indicating that the effects of A are 
dominant in the future late stages of the evolution of A-LTB configurations in which 
the worldlines of expanding dust layers are infinitely inextensible. As a consequence, 
with the exception of the unstable loitering models, all A-LTB models endowed with 
a late time evolution are compatible with the cosmic no hair" conjecture [21 [25], with 
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the asymptotic de Sitter state identified with the critical point FA. For dust layers 
collapsing from infinity (L — ^ oo), the phase space trajectories emerge from a past 
attractor which exactly coincides with the point FA. 



Seen from the inhomogeneous projection (bib) (panels (b)), the evolution 



towards FA of phase space trajectories of perpetually expanding layers contains extra 



information not available in the projection (51a|. Trajectories of layers with negative 
spatial curvature (figures [2j|3]) approach the Minkowskian saddle SP2, whereas layers 
with positive curvature approach the saddle SP4 (figures |4j [s] [s] and [9]) associated 
with spatially flat conditions. The approach to SP2 clearly indicates an intermediate 
low density state, and also occurs in trajectories of the FLRW phase space (notice 



that SP2 is projected by (51a) into the single saddle of this phase space). However, 
the approach to SP4 has no equivalent in trajectories of the FLRW phase space of 
layers evolving towards the future attractor. 

In configurations containing dust layers that bounce and collapse {zq changes sign: 



figures l5||9|) we used the system ( lla)-( lid) to obtain the full phase space trajectories 



In all cases the phase space trajectories in the homogeneous projection start in the 
past attractor FA, reach infinite values tl^'^^ — )■ oo as = and return along the 
same trajectories to the same attractor (which is now a future attractor). However, in 
the inhomogeneous projection (panels (b) of figures [5]and|6]) the curves do not return 
to the attractor along the same trajectories: if the pattern is expansion/collapse (figure 
[H]) they reach 6^^'^ oo as Zq — and return to the attractor from (5*^^^ — ?> — cxi. If the 
pattern is collapse/bounce (figure [i]), the trajectories reach (5^^^ — > — oo and return to 
the attractor from S^^^ — >■ oo. The same effect would happen in the mixed pattern 
configurations displayed in figures [8] and [9j but was not included in the plots as the 
full curves would make a very messy pattern. This difference in the behavior of curves 
emerging/returning from/to an attractor is an effect inherent in the inhomogeneity of 
LTB solutions. 

The bouncing and loitering models deserve a separate mention, as they have no 
equivalence in the case A = 0. In the bouncing models (figure |6|, layers initially 
collapse from infinity, hence their past attractor coincides with the future attractor of 
perpetually expanding configurations (figures [2]^. In the loitering configuration that 
we examined (figure [?]), the phase space trajectories emerge from the past attractor 
PA, but as the layers become asymptotically static — >• the trajectories evolve 
towards d^'^\Cl^^\S^^^ — >■ oo. In the homogeneous projection (panel (a) of figure [?]) 
the trajectories are indistinguishable from those of figure [5] but in the inhomogeneous 
projection (panel (b) of figure [?]) the curves reach (5'^^ — >■ oo along different paths and 
converge towards (5^^) 0. This effect is an inherent feature of inhomogeneity and 
does not occur in FLRW phase space diagrams for loitering models. It is worthwhile 
remarking that loitering models do not evolve towards an asymptotic de Sitter state 
associated with the future attractor FA. 

The dynamics of A-LTB models is not much different (qualitatively speaking) 
from that of the homogeneous FLRW models that they generalize. Essentially, the 
past attractor is a state close to spatially flat FLRW and self-similarity, while the 
future attractor is a de Sitter state dominated by A (in agreement with the "cosmic no 
hair" conjecture). The two most significant effects that arise from the inhomogeneity 
of A-LTB models are: 

(i) the possibility of accommodating wholly different kinematic patterns in the same 
spacetime configuration (see figures ^nOl). However, even when the kinematic 
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patterns of all layers are qualitatively similar, each unique A-LTB spacetime can 
be understood as an inhomogeneous model made up of some sort of superposed 
FLRW models (as the phase space trajectory of each comoving worldline can be 
projected by (51a| into the phase space trajectory of a unique FLRW model), 
(ii) the asymmetry in the evolutions from a past attractor and a bouncing point and 
from that point to the same attractor, which becomes a future attractor (see 
panels (b) of figures [s] and [6| 

However, in spite of the simplification involved in assuming spherical symmetry and 
geodesic motion, the present dynamical system study clearly shows that even this 
idealized level of inhomogeneity does provide extra degrees of freedom that can be 
be very handy for constructing useful toy models of astrophysical and cosmological 
inhomogeneities. 
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